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Motivated by recent experiments, we develop a model for DNA toroids under external tension. 
We find that tori are the equilibrium states for our model up to a critical tension, above which they 
become only metstable. Above this tension, we find a cascade of transitions between discrete toroid 
states that successively lowers the winding number, until the ground state (rod) is reached. In this 
process, this model predicts a nearly constant force plateau as a function of extension, in agreement 
with experiment. 



I. INTRODUCTION 

It has long been recognized that the conformation of 
polymer chains depends on the solvent properties of the 
environment [U, 0]- In particular, polymers in poor sol- 
vent conditions effectively attract each other in an at- 
tempt to exclude the solvent, forming collapsed struc- 
tures that minimize surface contact with the solvent. 
For flexible polymers, this leads to compact globules of 
roughly spherical shape, whose kinetic pathway has been 
shown to involve the formation of a pearl necklace and 
gradual diffusion of large pearls to the chain end [1, 0, Q ■ 

In the case of semiflexible polymers such as DNA, 
which exhibit a substantial bending stiffness, the ener- 
getic penalty for bending causes spherical globules to 
be energetically disfavored. The apparent equilibrium 
states for these polymers have been shown to be rings 
or toroids ^Q], as these structures balance the tendency 
for the polymer to condense due to effective polymer- 
polymer attraction with the tendency to minimize cur- 
vature due to bending stiffness. These condensed states 
have been studied theoretically EiPil p[-li^-lfl1 - fl2l . iT3t . ob- 
served in experiments [l4 j_ m [iTrilSl. Il9l . and shown 
by computer simulationlioliinilliinij^ 

In an effort to understand the dynamics of toroid for- 
mation, recent experiments have explored the conden- 
sation of DNA under tension [H, HI, [13, [H, 13. Mo- 
tivated by these experiments we analyze theoretically a 
hierarchy of tori states and explore their equilibrium and 
metastable structure under tension, as well as transitions 
between toroid states. We find a sequence of metastable 
tori under tension. Furthermore, we find that for winding 
numbers larger than approximately 10, a nearly constant 
force plateau emerges, which agrees well with recent ob- 
servations, as illustrated in Fig. [T] 

In Sec. m] of this paper, we first define a simple, non- 



thermal model that incorporates the essential physical ef- 
fects believed to give rise to DNA toroids: (1) the bend- 
ing rigidity and (2) the effective attractive interactions 
between DNA segments, such as can arise in the pres- 
ence of multivalent ions Q. We summarize the relevant 
experiments in Sec. IIIII We then study the equilibrium 
and metastable states of this model, as well as transi- 
tions among these states in Sec. IIVI We conclude with 
a discussion of the implications of this model, and the 
relationship of our results to the experiments. 
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FIG. 1: Measured force-extension curve on condensed DNA. 
A lambda-DNA molecule (48502 base pairs, or ~ 16.4/im in 
length) was stretched using optical tweezers in the presence 
of spermine4-|- (black line, averaged to 8 Hz, raw data at 128 
Hz shown in light gray). In contrast to the force-extension 
curve of uncondensed DNA (solid gray line) we find a force 
plateau of approximately 4 pN that persists throughout the 
curve. In this plateau small steps can be discerned, signifying 
unwinding of loops from a toroidal DNA condensate. 



II. THE MODEL 
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The first step in examining the equilibrium and 
metastable structure of semiflexible polymer condensates 
under tension is to identify and calculate their energy. 
We model the conformational energy of a toroid as in 
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Schnuur, et al. [l^l, where we assume integer wind- 
ing number toroids with a single radius of curvature at 
zero temperature. Our calculations describe a simplified 
model of tightly packed filaments of vanishing thickness. 
We do not take into account any winding defects due to 
topological constraints or variations in curvature due to 
filament thickness. 

We write the Hamiltonian of our system as a sum of 
bending and interaction terms 

H = Hbend + Hint, (1) 

where the bending term models the energy of the curva- 
ture of the major radius of the torus and the interaction 
term models the self-attraction of the polymer, or equiv- 
alently its poor solvent environment. The bending term 
can be straightforwardly calculated, as the energy of the 
idealized chain described is simply that of a series of cir- 
cular rings, given by 

where C(s) is the curvature, L is the total filament length 
in the torus, R is the torus radius, N is the torus winding 
number, and k is a bending stiffness constant. We note 
that since we consider only integer winding number tori, 
L = 27rNR, so we can write ([2|) in terms of a single 
quantity, L. 

For the interaction term, we assume a dense structure, 
in which filaments pack tightly in their plane perpendic- 
ular to their local axis. This suggests a simple, hexag- 
onal packing of the filaments. With such tight packing, 
we assume that the interactions are only of the near- 
est neighbor type. In this limit, the filament can be 
thought of as having six possible binding sites per unit 
length, which can either form a DNA-DNA bond with 
another section of the filament or can be exposed to sol- 
vent. Bundling occurs when the attractive interactions 
are sufficiently strong. In order to calculate the inter- 
action term, we define the number of bonds per cross 
section to be nt = ^ where Us is the number of 

solvent-exposed sites. We divide by two to avoid dou- 
ble counting, as a DNA-DNA bond is equivalent to the 
merging of two binding sites on neighboring filaments. 
Our interaction term can thus be written as 

Hint = -7 / ds nb{s) (3) 
Jo 

= -3jL + - ds ns{s) 

^ Jo 

= -aaL + "f—L, 

where 7 is a surface tension parameter that characterizes 
the energetic cost of of solvent-exposed DNA and ajv is 
the so-called coordination number , which is equal to 
half the number of solvent exposed sites per unit length 
along the torus circumference. We use the coordination 
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TABLE I: Coordination numbers for tori winding numbers 
1-24, with magic numbers in boldface. 



number to enumerate these sites, and it can be found by 
subtracting the number of filament-filament bonds from 
3A^. As an example of this scheme, consider the cases 
N = 5 and 10: for five filaments there are seven bonds, 
resulting in a coordination number of 8 (see Fig. [2]), while 
for ten filaments there are 19 bonds, resulting in a coordi- 
nation number of 11. We then multiply the coordination 
number by the interaction parameter, 7, to obtain the 
surface energy per unit length of a bundle. Table [J lists 
the first 24 coordination numbers. We replace the 37L 
term in the second line of Eq. [4]with jLai to emphasize 
the physical meaning of this term. The jaiL term comes 
from the difference in surface energy between N strands 
of unbundled filament and N bundled strands, and it 
reflects the physical tendency of the torus to minimize 
unsatisfied bonds through bundling. 

The packing of filaments in a hexagonal crystal brings 
about particularly stable toroids for certain winding 
numbers, as noted previously in Refs. [13, [3 ■ This sta- 
bility can be attributed to a high degree of hexagonal 
symmetry and the resulting low surface energy. From Ta- 
ble [J we can see that the difference between subsequent 
coordination numbers is either or 1 (for A^ > 2). 
For instance, in the case of the 5-torus bundle shown on 
the right of Fig. [21 the addition of a sixth filament to the 
bundle can satisfy no more than two bonds, resulting in 
no fewer than four additional unsatisfied bonds and an 
increase in the coordination number by one. If this sixth 
filament is added just above the filament on the right, 
then the addition of a seventh filament directly above the 
center results in no increase in the coordination number, 
since three bonds can be satisfied. Here, the result is a 




FIG. 2: Sketch of cross section of A'^ = 5 torus, with hexagonal 
filament cross section. The perpendicular lines represent bond 
between two sites on neighboring filaments. 
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symmetric, compact cross section that we refer to as a 
filled shell. In such cases, where = ct /y i , we refer 
to A'^ as a magic number, following Refs. |12l . |l3| . Fig- 
ure [3] shows toroid cross sections for the first seven magic 
numbers. The magic numbers, up to iV = 24, are as fol- 
lows: N ^ 7,10,12,14,16,18,19,21,23,24. If, instead, 
QfAT = ajsr-i + 1, then the smaller — 1 torus is favored 
by both the interaction energy, as well as the bending 
energy. 

As winding number goes up, we have a higher density 
of magic numbers. This can be understood by the in- 
crease in edge vs corner filaments in the filled shells for 
large N. For instance, the = 16 structure in Fig. [3] 
can be obtained by removing three edge filaments from 
the = 19 structure. All but the last one of these fil- 
aments satisfy three bonds, corresponding to no change 
in a, while the final corner filament satisfies only two 
bonds, corresponding to a reduction in a. Thus, we find 
sequences of increasing length of successive magic num- 
bers, although equilibrium tori are only found for the 
largest A^ in each sequence. 

These sequences can be identified as follows. Certain 
winding numbers correspond to perfect hexagons, such as 
N = 7 and N = 19 in Fig. [21 which we call supermagic 
numbers, as in Ref. [ll] . Since these hexagons consist of 6 
equilateral triangles of length fc = 1, 2, 3, . . . filaments on 
each side plus one filament in the middle of the hexagon, 
the supermagic sequence is given by A'^ = 3fc(fc + 1) + 1. 
The coordination number in this case is ajv — 3(2fc -I- 1). 
Thus, the difference between successive supermagic a is 
6, although the difference in A" is 6fc, where k corresponds 
to the larger A^. Hence, the length of each sequence of 
N with the same coordination number is fc, on average. 
In fact, as suggested by Table 1, the actual lengths of 
the sequences of successive A^ with equal coordination 
number are given by fc — 1, fc, k, fc, k, k + 1. 

As noted before, semiflexible polymer condensates 
generically form rings or toroids from the competition 
between their tendency to minimize surface area due 
to short-range attractive forces and their tendency to 
straighten out due to their substantial bending stiffness. 
Balancing these two effects, i.e., setting n/L ^ -fL, lets 
us define a natural length scale for our problem, which 
we call the condensation length, Lc — \/ k/^. Physically, 
this length is the approximate length scale at which we 
expect condensation to occur. Below this length DNA 
will rarely self-intersect and thus rarely condense, while 
above it a DNA filament will self-intersect many times 
and thus form collapsed, intermediate structures. 
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FIG. 3: Sketch of toroid cross sections for the first seven magic 
numbers. Note that while all have high degrees of symmetry, 
the N = 7 and A'' = 19 cross sections correspond to perfect 
hexagons, and we dub these winding numbers supermagic. 



We can also define an analogous energy scale, the con- 
densation energy, Uc — ^J^. Given these scales, we 
can present our conformational energies in dimension- 
less units, with physical values of length and energies 
normalized by their condensation values: F-^ = Un/Uc, 
where Un is the conformational energy of an A^ torus, 
and A = L/Lc. The presentation of our results in di- 
mensionless units clarifies the relevant parameters in our 
theory, namely the stiffness constant k and the interac- 
tion parameter 7. 

Combining our expressions for bending and surface en- 
ergy and normalizing by the condensation energy gives 
us the dimcnsionlcss free energy for an A'^ torus 

Fn = + - ai) (4) 

A. Equilibrium Torus States 

Plotting Eq. ^ for various A^ gives a family of curves 
which tend to negative infinity. While there is no definite 
global minimum in the free energy for all lengths, for a 
specific reduced length A, there is an associated optimal 
winding number N(X), which corresponds to the equilib- 
rium state at that length. For A < 47r, the rod, A^(A) — 0, 
is the equilibrium state. From A^ = 2 to A^ = 7, every 
state is an equilibrium solution except for A^ = 6, which 
we expect from our discussion of magic numbers. (These 
results differ somewhat from Ref. [13], since we focus on 
only tori of integer winding number.) 

Above N — 7, all equilibrium states are magic number 
states, though not all magic number states are equilib- 
rium states. Instead, only the largest in each sequence 
of consecutive magic numbers winding numbers corre- 
spond to equilibrium tori in the absence of tension. This 
can be understood as follows. For A^ > 12, the ratio 
aj^/N, and therefore the interaction energy, has local 
minima at each value of A^ such that both a at = Q!Ar_i 
and ayv+i = oln + 1- This forms a sequence of winding 
numbers A^ = 14, 16, 19, 24 . . . for which particularly sta- 
ble tori are expected. As can be seen in Tab. |T1 these 
all correspond to magic numbers. Including the effect 
of bending energy, which always favors smaller winding 
numbers, consistent with Ref. p^ . we find stable tori 
for this sequence, as well as for the other magic numbers 
A^ = 7,10, 12. 

It has been noted before [H, [3 ^^at the equilibrium 
radii of the tori do not increase monotonically as a func- 
tion of reduced length. In fact, as reduced length is in- 
creased, the radii of subsequent equilibrium tori states 
are marked by discontinuous jumps. Fig. |3] shows the 
reduced radius, p, of the equilibrium states as a func- 
tion of reduced length, up to A^ = 19. The discrete 
transitions between the radii of different winding num- 
ber is again an effect of the hexagonal packing, which 
creates islands of stability for certain winding numbers. 
An equilibrium toroid grows in radius until it reaches 
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a contour length at which the next equiHbrium winding 
number is favored, at which point it transitions to this 
state. The extra length needed for the additional loops 
of the higher winding number torus drives the toroid to 
take on a smaller radius for the same contour length. For 
comparison, we also show as the dashed line in Fig. |4]the 
prediction based on a continuous approximation valid for 
large N [9, 13]. 




FIG. 4: Plot of reduced radius vs reduced length for equilib- 
rium tori up to A'^ = 19. The dashed line indicates the pre- 
diction based on an approximation valid for large A'^ [l^ . 



B. Tori Under Tension 

The dimensionless free energy expression given in 
Eq. (IH) is a measure of the internal energy of the toroid. 
At zero temperature this is analogous to the Helmholtz 
free energy of our system. This free energy depends on 
the reduced length of the filament, or the volume of our 
essentially ID system. Since we also want to consider 
the effect of tension, which as a force variable plays the 
role of pressure in a classical thermodynamics analogy, 
we perform a Legendre transform on the internal energy 
of the toroid to get our energy expression in terms of 
force 
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(5) 



where we equate ^^^jf- = /■ Here, Eff is analogous 
to the Gibbs free energy, where tension is the control 
variable. 

The dimensionless form of the Gibbs free energy is ob- 
tained again by normalizing Eq. ([5]) by Uc 
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where t = f is the dimensionless tension. Equation ([6]) 
represents the dimensionless result for the toroid's con- 
formational energy under tension that we will generally 
be referring to when we discuss toroidal energy. 



III. REVIEW OF EXPERIMENTS 

Over the past few years, a number of single- molecule 
experiments have probed the mechanics of DNA conden- 
sation under tension [2^, [2^, |23, Ull, [s^]. Generally, in 
these experiments a single DNA molecule is stretched 
and relaxed using optical or magnetic tweezers. Under 
conditions appropriate for condensation, most of these 
studies reported a nearly constant force plateau of sev- 
eral pN for DNA extensions lower than ~85% of the full 
contour length. While exact numbers differ between ex- 
periments, similar qualitative behavior was observed for 
a wide range of condensing agents and concentrations. 

It has been generally believed that the force plateau 
regime consists of a continuous unraveling of DNA un- 
der tension. In order to test this hypothesis, wc mea- 
sured the force-extension relationship of a single con- 
densed DNA molecule with high resolution using opti- 
cal tweezers. Lambda-DNA (48502 bp) was attached on 
both ends to two optically trapped polystyrene beads [2^ 
and allowed to condense in the presence of 1 mM sper- 
mine4-|-. The condensed DNA was subsequently unrav- 
eled by displacing one of the beads. In these experiments, 
as shown in Fig. [U the force was allowed to vary freely, 
while the extension was varied by displacing the trap. We 
observed a roughly constant force plateau that was con- 
sistent with previous experiments. In detail, however, we 
found that this plateau consists of a sawtooth-like pat- 
tern, which suggests a step-like unraveling of the DNA 
under increasing extension. With the model above, we 
can account for these discrete steps in terms of jumps be- 
tween toroid states with different winding numbers. In 
addition, this model can also account for the nearly con- 
stant force plateau. 



IV. RESULTS OF THE MODEL 

We find, for finite filament length and zero tension, 
that tori states are, indeed, the equilibrium states within 
our model, as discussed above. The filament length dic- 
tates which torus state has the lowest free energy. For 
a given filament length A and winding number N, there 
exists a finite tension Tcrit at which the N torus begins 
to unravel. For tension r < (ai — on/N), the free en- 
ergy Gn strictly decreases with increasing A, resulting 
in a stable or metastable state in which the torus incor- 
porates the full polymer length. Even before this point, 
however, the free energy Gjv(''';A) may become greater 
than zero, indicating that the thermodynamically stable 
state is the extended polymer conformation. In this case, 
the tori are actually metastable states. 

With increasing tension r > (ai — aj^/N), a local 
minimum in Gn develops for A less than the full poly- 
mer length. Physically, this corresponds to a mechanical 
metastable state, or local energy minimum, in which a 
torus coexists with a segment of unwound, straight poly- 
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mer under tension. The condition for this to occur is 



where the derivative on the right is evaluated at the full 
polymer length. As the tension increases, more filament 
is pulled out of or unwound from the torus, which then 
shrinks in size. (Here, and throughout, we assume that 
the torus is able to relax by internal relative sliding of 
polymer.) As the tension increases and the torus shrinks 
in size, the increased bending energy eventually results in 
destabilization of the A^-torus relative to tori of smaller 
winding number. We identify below a series of transitions 
under tension to tori of smaller winding number. 

A. Transitions 

As tension is increased from zero, the slope of the free 
energy curve for each N increases. With increasing ten- 
sion, the N[X) torus initially in equilibrium in the ab- 
sence of tension for a given total length A remains the 
lowest energy state until a critical tension is reached. At 
this critical tension, Tcru, the asymptotic slope of the free 
energy becomes zero for large A. For tensions above this 
critical tension, the N{X) torus becomes only metastable, 
as the (TV = 0) rod is now energetically more favorable 
(with a free energy of zero) , and is thus the equilibrium 
state of the system. As the tension continues to increase, 
a local minimum of the free energy develops and starts 
to shift to lower values of A. Once this minimum shifts to 
values of A less than the full length of the DNA strand, 
then the metastable state consists of a compact N torus, 
with a segment of filament pulled out of the torus — i.e., 
the torus begins to unravel. With increasing tension, as 
more of the filament is pulled out of the torus, larger N 
tori become unstable to tori with smaller winding num- 
bers. This unraveling process is sequential, with transi- 
tions to smaller and smaller values of N as the tension 
increases. 

With an eye toward addressing the experiments in 
Refs. [H,|ll|, we consider a process in which the reduced 
length A in the torus is controlled and slowly reduced, 
while the tension is allowed to vary. As the N torus is 
slowly unravelled, the tension r will increase. However, 
as T increases, we expect at some point to develop a local 
minimum of Eq. ^ ioi N — 1 that becomes less than or 
equal to that of the (metastable) N toroid. Once a tran- 
sition to the — 1 toroid occurs, if the length A is fixed, 
then the tension will fall as a new (metastable) A^ — 1 
toroid is formed at A. This describes most of the unravel- 
ling transitions, at least for large A, where a^-i — olm- 
However, given the discrete nature of the coordination 
number oat, it can happen that the A — 1 state is itself 
unstable to the A — 2 state at A. This occurs for A > 9 
whenever aN-2 < aAf-i, i-e., when A — 2 is a magic 
number. Thus, the A — 1 state can be expected to make 
a transition to the A — 2 state: in the limit of a slow 



unravelling of the torus, the A — 1 state is skipped. The 
resulting sequence of states and corresponding tensions 
vs reduced polymer length A is shown in Fig. O Note that 
the extension of polymer pulled out of the torus varies 
inversely in A, meaning that metastable branch corre- 
sponds to a stable force-extension relation, in which force 
increases as more polymer is pulled out of the torus. We 
find, interestingly, for A greater than about 7, a nearly 
constant force plateau at t w 2.5. This is consistent 
with several recent experimental observations of a nearly 
constant force for tori under tension [2^ [H, [H, [2§| , as 
illustrated in Fig. [1] 

The force plateau can be interpreted as the average 
force needed to pull off a loop from the torus. Since 
the number of bonds broken when pulling off a loop is 
either 2 or 3, depending on whether or not the toroid 
has a magic winding number, it's not surprising that we 
should see a force plateau at r 2.5. As we go to higher 
and higher winding number, however, the force plateau 
asymptotes to 3. We can see why this is by considering 
the average number of bonds per filament as a function 
of winding number. In Fig. [SI the most weakly bound 
filaments are the corner ones, which satisfy only three 
bonds. The remaining filaments along one edge adjacent 
to this corner can also be removed at the cost of just 
three bonds each, until the final corner filament along 
that edge is reached, the removal of which involves the 
breaking of just two bonds to form the magic number 
bundle with the next lowest coordination number (A = 
16 and A = 33 for the bundles in Fig. [6]). With increasing 
winding number, the fraction of filaments forming three 
bonds increases, and r — > 3, although this convergence is 
slow. 
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FIG. 5: Tension vs reduced length for preferred tori states 
from A = 2 to A = 24, only plotted at tensions and lengths 
where we expect them. We note that only magic winding 
numbers appear for N > 7 and that a rough force plateau 
appears at r « 2.5. Transitions between winding numbers are 
characterized by discontinuous jumps in tension and length. 
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FIG. 6: Cross sections of JV = 19 and iV = 37 tori with the 
hghter circles corresponding to corner filaments. As winding 
number increases the number of corner filaments always re- 
mains 6, so the overall density of corner filaments goes down 
for large tori. 



B. Large N Behavior 
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FIG. 7: The range of reduced tensions r predicted for the 
metastable states, beginning with the 4^3 transition, using 
the large- A approximation of Eq. (|9} . 



For large winding number N, the number of exposed, 
unsatisfied bonds at the perimeter of a torus cross section 
increases as Vn, so that aN Vn. The prefactor here 
is easy to calculate for perfect hexagons, as illustrated in 
Fig. [6l As noted in Sec. |TT1 these occur for N = 3k{k + 
1) -I- 1, where k = 1,2,3,..., for which the coordination 
number is 

aN = 3(2fc + 1) = V3 {4:N - 1). (8) 

This actually represents a lower bound on a^Vi in gen- 
eral, since less symmetric cross sections have increased 
surface-to- volume or circumference-to-area ratios. As we 
are interested in the large N behavior, we will approxi- 
mate 

ttAT ~ 2V3N, (9) 

as in Ref. [Tsj . 

We use this large N approximation to determine the 
metastable states and transitions between them, as we 
have done in the previous section. Specifically, we con- 
sider fixed but decreasing condensed length A. In Figs. [7] 
andm we indicate the predicted sequence of metastable 
states and corresponding values of tension r and toroid 
size p in reduced units. We find many of the same qual- 
itative features in this large N approximation as were 
found in the previous section. In particular, we find an 
apparent force plateau, much as in Fig.O Although both 
models predict the same asymptotic convergence of the 
metastable tension r — > 3 for large A, where A^ is also 
large, we see that the approximate value of from 
Eq. ([9]) yields a consistently smaller value of the tension 
than for the discrete model in Sec. IIV Al This can be 
understood as follows. In the discrete model, we account 
for the hexagonal packing of filaments, which results in a 
sequence oi aN in which there are discrete jumps in ajv 
(at non-magic numbers), between which ranges of con- 
stant aN are found. For transitions from A^ to A^ — 1 
toroid states, the A^ — 1 state is thus destabilized in the 
discrete model as compared to the large N approxima- 
tion, when aN = aN-i- This enhanced relative stability 



of the A'' toroid means that the tension r is larger at 
the transition. For the same reasons, the corresponding 
toroid sizes are smaller in the discrete model than in the 
large A^ model. Nevertheless, the general features, and 
especially the force plateau, are seen for both models. 



C. Energy Barriers 

We estimate the energy barrier between A^ and N — M 
tori {M < N) first by calculating the energy difference 
between an A^ torus and N — AI + M tori (see Fig. 
This is not to suggest that Fig. [H] represents the actual 
reaction pathway for the transition between tori: calcu- 
lating the energy in this way, by considering the separa- 
tion of full loops from the original torus, can only pro- 
vide an upper bound on the transition energy. From our 
model, and as suggested by experiment, we expect that 
toroids make step-wise jumps between different winding 
numbers. 
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FIG. 8: The range of reduced torus radii p predicted for 
metastable states, beginning with the 4 ~+ 3 transition, using 
the large- A approximation of Eq. ® . 
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FIG. 9; Sketch of TV, iV-M and M tori. Physically we expect 
the energy barrier between these to be dependent only on 
surface energy difference, as the bending energy and tension 
are unchanged. 



Intuitively, we expect the transition to depend primar- 
ily on the difference in surface area exposed to solvent, 
and not on the bending energy, since the torus radii vary 
only weakly with N. In what follows, we shall assume 
that the radii are constant. We can see from Fig. [9] that 
the following relations hold, since the total amount of re- 
duced length, Aat, is unchanging and all the loop radii 
are equal: Aat — Xn-m + Am, where \n-m = ^7?^^ A 



N 



N 



and Am = ^Xn- Here, Aat is the length in the N torus, 
Xn-m the length in the N — M torus, and Am the length 
in the M torus. 

With these identifications we can straightforwardly 
calculate the energy difference between the N torus and 
the N — M + M tori, from the energy expression given 
by Eq. ® 



AE 



N,N~M 



+tXn-m + 



Ajv-M 



+ Xn-m{ 



aN-M 
N ~ M 



Oil) 



-^ ^ ^Af(-r7- - cti) + tXm 

Am jW 



Ajv 
Aat 



N 



(10) 



which depends only on the surface parameters, as ex- 
pected. We calculate Xn in the above equation as follows; 
starting from an N torus we assume downward sequential 
transitions in winding number, as previously described, 
occurring for Xn where the metastable N and A^— 1 states 
have equal free energy. For unstable tori, we assume a di- 
rect transition to the subsequent winding number at the 
same value of A. It is not our purpose to calculate realistic 
rates of such transitions. Instead, we focus on identifying 
the most relevant sequence of metastable torus states in a 
fixed-A ensemble, corresponding to experiments in which 
the total length is the control variable. Figure fTOl shows 
the energy barriers between N ~* N ~1 transitions, from 
iV = 24 down. We find that transitions from N to N — M 
for M > 1 have substantially higher energy barriers, and 
are thereby strongly suppressed. 

As noted above, it is unlikely that the transitions 



sketched in Fig. [9] represent the real reaction pathway 
between torus states. Nevertheless, our estimates sug- 
gest that transitions corresponding to a change in wind- 
ing number by more than 1 are strongly suppressed. For 
AA^ = 1, however, it is also possible that DNA peels 
off continuously from the torus. This will, in general, 
induce additional bending as the DNA is pulled off the 
torus. (We assume twist is relaxed.) In fact, a bend 
of approximately 90 degrees is expected on simple me- 
chanical grounds, since the force in the plateau region 
in Fig. is comparable to the total binding energy per 
unit length of filament, as discussed in Sec. IIV Al Un- 
der a force /, the radius of curvature r for such a bend 
can be estimated by balancing the total bending energy 
7rK/(4r) and the virtual work against the applied tension 
(7r/2 — l)r/. This yields an optimal radius of curvature 
r = ^TTK/[f{2Tr — 4)], and a total energy for two such 
bends of •\/T7r(27r — 4) « 4.2 — 4.6 in reduced units for 
r « 2.5 — 3. This is about a factor of 2-3 smaller than 
the estimates from Eq. (jlOp . 



Other effects may change the energy barriers as well, 
such as next-nearest neighbor interactions. As these in- 
teractions are effectively attractive, they will favor more 
compact structures, thus lowering the length needed for 
an N torus (i.e. shifting minima in the energy to lower 
values of A). This will lower the energy barriers given 
above for nearest neighbor interactions, but there will be 
an additional contribution to the energy barrier from the 
attraction of next-nearest neighbors that must be over- 
come to pull off a loop. Without explicit inclusion of 
these effects in our model it is unclear what net effect 
these interactions will have on energy barriers. 



40 - 



30 - 



N 



FIG. 10: Estimated energy barriers (dimensionless) vs torus 
winding number, using Eq. (|10[l . for all A'^ to A'^ — 1 transitions 
between N = 2 and N = 24. The energy barriers are between 
TV and A^ — 1 tori, e.g., the energy barrier plotted at TV = 12 
is the energy barrier between TV = 12 and TV = 11. 
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V. DISCUSSION 

Of the two parameters in our model, the bending 
stiffness k is known to be approximately k = kTip ~ 
SO/cT-nm, where ip is the persistence length of DNA. 
The interaction parameter 7 is expected to depend on the 
counter ions present in solution, and it is unfortunately 
not known. In order to estimate 7, we use the experimen- 
tally observed plateau value of the force in Fig. [TJ This is 
consistent with the plateau we find in Fig. [Sj By match- 
ing the measured force plateau with our (dimensionless) 
tension r — 7/7 we estimate 7 ~ 1.6pN. This allows us 
to estimate the condensation length, Lc — llnm, and 
condensation energy, Uc — 4:.5kT. This corresponds to 
energy barriers of order 20fcr, using the lower estimates 
at the end of Sec. IIV CI for AiV = 1 transitions. By con- 
trast, we find substantially larger barriers of more than 
50fcT for transitions as sketched in Fig. [9] for AA^ > 1. 

Figure [5] shows a typical transition length (i.e., the dif- 
ference in reduced length between an N torus and an 
TV — 1 torus at constant tension) of ~ 10, suggesting 
typical torus loo p si zes of ^ 112nm, a value larger than 
reported in Ref. [23] by a factor of about 2-3. Given the 
simplicity of our model, being off by such a factor is per- 
haps not so bad, as there are many effects that we have 
not taken into account. One notable effect that we've ne- 
glected is next-nearest neighbor interactions. Such inter- 
actions will have the same attractive tendency as nearest 
neighbor bonds, and thus will tend to favor more com- 
pact structures, translating into smaller loop sizes. 

Another possible effect is stability of the toroids be- 
yond the transition regime. We have predicted the loop 
sizes from our model assuming immediate transitions 
once in the transition regime; in reality toroids may 
not transition immediately. At higher tensions there are 
smaller differences in reduced length between subsequent 
toroid force-extension curves. If, for example, we instead 
assume a force plateau at the higher-than-expected value 
T — 3 (which raises Lc to 12nm), and average over the 
same range of winding numbers, we get an average step 
size of A ~ 6.5, translating into a loop size of 78nm. 

One very important effect that must be taken into ac- 
count in comparing with the measured extensions in the 
experiments is the finite extensibility of the uncondensed 
DNA strand, which is expected to be well approximated 
by a worm-like chain [28', '291, [s^ . We can account for this 
by considering a thermally fluctuating filament in series 
with the torus. The length of this filament is equal to 
the full contour length Lq of the DNA in the absence of 
condensation, minus the length L contained in the torus. 
In dimensionless form, the extension ^ = x/Lc of the free 
DNA in series with the toroid is approximately given by 



where the full contour length of uncondensed DNA is 
given by ALc, of which XLc is contained in the toroid. 
In Fig. [TT] we show the combined force-extension curves 



4 



3 



T 2 



























N=2^ 






1000 


1100 


1200 1300 



FIG. 11: The predicted force (in units of 7) as a function 
of the apparent extension ^ (in units of Lc) of DNA toroids 
in series with freely fluctuating uncondensed DNA. Here, we 
have used 16.4/im as the full contour length and Lc = llnm. 
A few of the individual force-extension curves are labeled by 
the corresponding winding numbers < 10. The thin con- 
tinuous curve represents the force-extension curve of the bare 
DNA. 



of the the various toroid states in Fig. [5l corresponding 
to Vc — 4:.5kT and the full contour length of 16.4^m in 
Refs. [H, [2^. Interestingly, this model predicts mul- 
tiple metastable states with winding numbers in the 
range of approximately 2-10 near the transition to the 
fully unraveled toroid. This may explain the signifi- 
cant hysteresis reported for forward (extension) and re- 
verse force-extension measurements near the transition 
between toroid and fully extended DNA [2^, [2^. As 
the fully extended state is allowed to condense when re- 
ducing the extension, the force drops and toroid states 
with small TV are expected to form at a smaller extension 
than required, e.g., for the transition from ~ 7 to fully 
extended DNA. 

The simple model developed here is able to capture a 
number of features observed in the experiments. We find 
theoretically that toroid unraveling under tension occurs 
via a series of discrete transitions, as observed in Refs. 
[H, [2§|. In addition, this model provides an explanation 
for the fact the approximately constant force plateau for 
DNA condensations under tension, as reported in a num- 
ber of experiments [H [H, [13, [H, [ll] . 
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